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The production of isolated photons in deep-inelastic scattering ep — > ejX is measured 
with the HI detector at HERA. The measurement is performed in the kinematic range of 
t^^ , negative four-momentum transfer squared 4 < Q^ < 150 GeV^ and a mass of the hadronic 

^-^ ' system Wx > 50 GeV. The analysis is based on a total integrated luminosity of 227 pb~^. 

The production cross section of isolated photons with a transverse energy in the range 
^ . 3 < Ej, < 10 GeV and pseudorapidity range —1.2 < 7]'~' < 1.8 is measured as a function 

^ [ of Ej,, if I and Q^. Isolated photon cross sections are also measured for events with no jets 

or at least one hadronic jet. The measurements are compared with predictions from Monte 
Carlo generators modelling the photon radiation from the quaik and the electron lines, as 
well as with calculations at leading and next to leading order in the strong coupling. The 
predictions significantly underestimate the measured cross sections. 
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1 Introduction 

Isolated photons originating from the hard interaction in high energy collisions involving hadrons 
are a sensitive probe of perturbative Quantum Chromodynamics (QCD) [1,2], as the photons 
are largely insensitive to the effects of hadronisation. A good understanding of the Standard 
Model (SM) production mechanism of isolated photons is also important for searches of new 
particles decaying to photons at hadron colliders. 

The production of isolated photon^Jhas been studied at various experiments. Cross sections 
measured in fixed target pN experiments (e.g. [3]) show a steeper decrease with photon trans- 
verse momentum, P^, than predicted by next-to-leading order (NLO) QCD calculations. The 
CDF [4] and DO [5] experiments at the Tevatron have measured the isolated photon production 
cross section in pp collisions. Whereas DO finds good agreement with a NLO QCD calculation, 
the CDF data show a somewhat steeper P^ dependence than predicted. Measurements of the 
photon production in e+e collisions have also been performed at LEP [6]. At HERA, prompt 
photon cross sections have been measured by the HI and ZEUS experiments [7-9] in photo- 
production, where the negative four-momentum transfer squared Q^ of the exchanged virtual 
photon is close to zero, and showed reasonable agreement with NLO calculations. An analy- 
sis of the isolated photon cross section in deep-inelastic scattering (DIS) with Q^ larger than 
35 GeV^ has been published by ZEUS [10]. 
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Figure 1: Leading order diagrams for isolated photon production in DIS. The upper diagrams 
illustrate isolated photon production by radiation from the electron line {LL), while the lower 
diagrams correspond to production via radiation from the quark {QQ, without the contribution 
from quark fragmentation). 



The measurement of isolated photons in DIS provides a test of perturbative QCD in a kine- 
matic range with two hard scales: the transverse energy of the emitted photon E]. and Q^. 



'Photons coupling to the interacting partons are often called "prompt" in contrast to photons from hadron 
decays or those emitted by leptons. 



Isolated photons in DIS are produced at lowest order (a^a°) as shown in figure [IJ Already at 
this order a jet can be produced in the hadronic final state in addition to jets associated with the 
proton remnant, due to the electron or photon recoil. 

The final state photon may be emitted by a quark {QQ) and by wide angle radiation from 
the lepton {LL). The interference contribution {LQ) is expected to be small. Since the photon 
and the scattered electron are well separated in the present analysis, low angle QED radiation 
is suppressed. The QQ contribution has two different origins: the direct radiation of a photon 
from the quark and the fragmentation of the quark into a jet containing a photon which carries a 
large fraction of the jet energy. This quark-to-photon fragmentation contribution is suppressed 
by the isolation requirement for the photon. 

This paper presents a measurement of isolated photon production in DIS e+j»^e + 7 + X. 
Photons are identified using a multivariate analysis of the shapes of the calorimeter energy 
deposits to reduce the background from neutral hadrons and their decay products. The photons 
are then used together with the other particles in the event, with the exception of the scattered 
electron, to reconstruct jets. The isolation of the photon is ensured by requiring that it carries 
at least 90% of the transverse momentum of the jet containing the photon. Isolated photons 
with transverse energy 3 < E'^ < 10 GeV and pseudorapidity —1.2 < if < 1.8 are selected 
in DIS events in the kinematic regime 4 < Q^ < 150 GeV^, inelasticity y > 0.05 and a mass 
of the hadronic system Wx > 50 GeV. The production of additional jets besides the photon 
jets in these events is also investigated. The current analysis significantly extends the kinematic 
range probed by the ZEUS measurement [10]. The results are compared to a recent leading 
order (LO), 0{a^a^), calculation [11, 12] and to predictions of the Monte Carlo (MC) models 
PYTHIA [13], simulating the QQ process, and RAPGAP [14] for the LL process. The cross 
sections for a photon plus at least one jet are further compared to a NLO calculation [15]. 



2 HI Detector 

A detailed description of the HI detector can be found in [16]. In the following, only detec- 
tor components relevant to this analysis are briefly discussed. The origin of the HI coordinate 
system is the nominal ep interaction point, with the direction of the proton beam defining the 
positive z-axis (forward direction). Transverse momenta are measured in the x-y plane. Polar 
(9) and azimuthal (0) angles are measured with respect to this reference system. The pseudora- 
pidity is defined to be 77 = — lntan(6'/2). 

In the central region (20° < 6* < 165°) the interaction region is surrounded by the central 
tracking system, which consists of a silicon vertex detector, drift chambers and multi-wire pro- 
portional chambers, all operated within a solenoidal magnetic field of 1.16 T. The trajectories 
of charged particles are measured in the central tracker with a transverse momentum resolution 
of (j{Pt)/ Pt — 0.005 Ft /GeV © 0.015. The forward tracking detector and the backward drift 
chamber (operated in 1999-2000) or backward proportional chamber (for 2003-2005) measure 
tracks of charged particles at smaller (7° < 6* < 25°) and larger (155° <6< 175°) polar angle than 
the central tracker, respectively. In each event the ep interaction vertex is reconstructed from 
the measured charged tracks. 



The liquid argon (LAr) sampling calorimeter [17] surrounds the tracking chambers. It has 
a polar angle coverage of 4° < ^ < 154° and full azimuthal acceptance. It consists of an inner 
electromagnetic section with lead absorbers and an outer hadronic section with steel absorbers. 
The calorimeter is divided into eight wheels along the beam axis, each of them segmented in 
into eight modules, separated by small regions of inactive material. The electromagnetic and 
the hadronic sections are highly segmented in the transverse and the longitudinal direction with 
about 44 000 cells in total. The granularity is larger in the electromagnetic part and increasing in 
both sections in the forward direction. For particles coming from the ep interaction region, the 
laterally projected cell size in the electromagnetic part varies between 5x5 cm^ in the forward 
and at most 7 x 13 cm^ in the central region. The longitudinal segmentation in the different 
wheels varies from three (central) to four (forward) layers in the electromagnetic and from four 
to six in the hadronic section. The first electromagnetic layer has a thickness of about 3 to 
6 radiation lengths for particles coming from the interaction region. Electromagnetic shower 



energies are measured with a precision of a{E)/E = \2%/\jE/ GeY © 1% and hadronic 
energies with a{E)/E = 50%/ ^/E]~GeV © 2%, as determined in test beam experiments [18, 
19]. In the backward region 153°<^<178°, particle energies are measured by a lead-scin- 
tillating fibre spaghetti calorimeter (SpaCal) [20]. 

The luminosity is determined from the rate of the Bethe-Heitler process ep — ^ ep'j, mea- 
sured using a photon detector located close to the beam pipe at z = —103 m, in the backward 
direction. 

DIS events at Q^ values up to 150 GeV^ are triggered by the energy deposition of the scat- 
tered electron in the SpaCal. For events with the scattered electron entering the SpaCal at low 
radii, additional trigger signals are required from the central drift chambers [21,22] and from 
the central proportional chambers [23-25]. The trigger efficiency for DIS events containing an 
electron in the Spacal angular acceptance with an energy above 10 GeV is greater than 98%. 



3 Monte Carlo Simulations 

Monte Carlo simulations are used to correct the data for detector acceptances, inefficiencies and 
migrations and to compare the measured cross sections with MC model predictions. 

The two generators PYTHIA [13] and RAPGAP [14] are used to generate events with pho- 
tons produced in the hard interaction. PYTHIA simulates the contribution of photons radiated 
from the struck quark (QQ). The contribution of photons radiated by the electron (LL) is sim- 
ulated using RAPGAP and denoted "RAPGAP rad." in the following. The small contribution 
from interference [11] is neglected. Both generators calculate the hard partonic interaction in 
LO QCD C(a^a°). Higher order QCD radiation is modelled using initial and final state parton 
showers in the leading log approximation [26]. The fragmentation into hadrons is simulated 
using the LUND string model [27] as implemented in JETSET [28]. The simulations use the 
CTEQ6L proton parton densities [29]. 

The measurements presented in this paper show that the data is well described by the two 
MC contributions if PYTHIA is scaled by a factor 2.3 and RAPGAP is not scaled. This com- 
bined "scaled signal MC" is used to correct the data, whereas the unsealed MC prediction 
("signal MC") is compared to the cross section measurements. 



As an alternative, the HERWIG [30] generator is used to model the QQ contribution. HER- 
WIG simulates the fragmentation into hadrons through the decay of colourless parton clusters 
and uses the equivalent-photon approximation for the incoming photon beam. Isolated photon 
production in DIS is derived approximately as Compton scattering between the photon and a 
quark. This approximation is not valid for Q^ above a few GeV^, therefore HERWIG is only 
used to estimate the systematic uncertainties due to the fragmentation model. 

The main SM background is due to photons produced in hadron decays in DIS events. It is 
modelled using the RAPGAP generator, with initial and final state radiation switched off. This 
contribution is denoted by "RAPGAP non-rad" in the following. 

The multivariate shower shape analysis used to identify the photons requires high statistics 
samples of shower simulations in the whole phase space of energy and pseudorapidity. Samples 
of events containing single particles are simulated. In each sample, corresponding to single 
photons or single neutral hadrons decaying to photons, the particles are uniformly generated in 
pseudorapidity and energy. These samples are generically named "single particle samples" in 
the following. 

All generated events are passed through a full GEANT [31] simulation of the HI detector 
and through the same reconstruction and analysis programs as used for the data. 



4 Event Selection 

The event sample used in this analysis was collected with the HI detector at HERA in the period 
1999 to 2005 at a centre-of-mass energy of 319 GeV. The corresponding integrated luminosity 
is 227 pb^^ In a first step, DIS events are selected with the scattered electrono measured in the 
SpaCal. In a second step, a subsample of DIS events with an isolated photon candidate in the 
LAr calorimeter is selected. 



4.1 Selection of DIS events 

DIS events are selected with the scattered electron identified in the SpaCal as a compact electro- 
magnetic cluster [32] with an energy Eg > 10 GeV and a polar angle Of. < 177°. Matching sig- 
nals in the backward tracking chambers are required for electron candidates with Ee < 18 GeV. 
The scattering angle of the electron is determined from the measured impact position in the 
backward tracking chamber, the position of the energy cluster in the SpaCal and the recon- 
structed primary vertex. 

Background from events at low Q"^, in which the electron escapes through the beam pipe 
and a hadron fakes the electron signature, is suppressed by the requirement that the differ- 
ence T.(E — pz) between the total energy and the longitudinal momentum be in the range 
35<S(_E' — pz) <70 GeV, where the sum includes all measured hadronic final state particles 
and the scattered electron. 



^The analysis uses data from periods when the beam lepton was either a positron or an electron. Unless 
otherwise stated, the term electron refers to both the electron and the positron. 



Non-ej9 background is removed by restricting the ^-coordinate of the event vertex to be 
within ±40 cm of the average vertex position and by requiring at least one good track in the 
central tracking system with the polar angle 30° < 9 < 150° and not associated to the electron. 

The energy Ee and polar angle 9e of the scattered electron candidate are used to reconstruct 
y and Q^ according toQ^ = 2 E^Ell + cos 6^) a.ndy = l- Ee{l - cos ^e)/(2^e ), where E^ is 
the electron beam energy. The events are selected in the kinematic region 4 < Q^ < 150 GeV^ 
and y > 0.05. 



4.2 Selection of isolated photon candidates and jets 

Photon candidates are identified as clusters in the electromagnetic section of the LAr calorimeter 
with a transverse energy 3 < E^ < 10 GeV and pseudorapidity —1.2 < r]'^ < 1.8 in the HI 
laboratory frame. The candidates are rejected if they are close to inactive regions between 
calorimeter modules [17] or if a track geometrically matches the electromagnetic cluster with 
a distance of closest approach to the cluster's barycentre of less than 20 cm. Neutral hadrons 
that decay into multiple photons, predominantly 7r° — > 77, constitute the main background. In 
most cases such decay photons are merged into one electromagnetic cluster, which tends to have 
a wider transverse distribution than that of a single photon. The transverse radiugj Rt of the 
photon candidate cluster is therefore required to be smaller than 6 cm. In addition, the invariant 
mass of the cluster, when combined with the closest neighbouring electromagnetic cluster with 
an energy above 80 MeV, must be larger than 300 MeV. This requirement rejects candidates 
that originate from 7r° decays with two photons reconstructed in separate clusters. Only events 
with exactly one photon candidate are accepted. Less than 1% of the events are rejected because 
more than one photon candidate is found. 

The mass of the final state hadronic system is reconstructed from the four-momenta of 
the incoming electron (pe) and proton (pp), the scattered electron (p'J and the photon candi- 
date (p^) as Wx = \/{Pe + Pp — P'e~ P-yY- The contribution from elastic Compton scattering 
(jo + e— >p + e + 7)is suppressed by requiring Wx > 50 GeV. 

Final state hadrons are reconstructed from deposits in the LAr calorimeter in combination 
with tracking information. Following the so-called "democratic" procedure [33,34], the photon 
candidate and the reconstructed hadrons in each event are combined into massless jets using 
the kx algorithm [35]. The algorithm is used with a P^-weighted recombination scheme and 
with the separation parameter Rq set to 1. Jets are selected with a transverse momentum of 
P|,^* > 2.5 GeV and a pseudorapidity in the range —2.0 < 77^"^* < 2.1. Due to the harder 
kinematical cuts for the photon candidate there is always a jet containing the photon candidate, 
called the photon-jet. All other jets are classified as hadronic jets. For hadronic jets the r]^^^- 
range is restricted to —1.0 < 77^"^* < 2.1. According to the MC simulation, the hadronic jets are 
well correlated to the partonic jets even at low transverse energies. To ensure isolation of the 
photon, the fraction z of the transverse energy of the photon-jet carried by the photon candidate 
has to be larger than 90%. This definition of the isolation requirement is stable against infrared 
divergences and thus well suited for comparisons with perturbative QCD calculations. The 
isolation requirement largely suppresses background from photons produced in hadron decays. 



^For a definition of the transverse radius see section |5] 
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The distributions of the transverse energy and of the polar angle of the isolated photon 
candidates are shown in figure [2] together with the MC predictions for the scaled signal and the 
background. The sum of the MCs describes the data well. 

The samples of events with either no hadronic jet or at least one hadronic jet are called 
"photon plus no-jets" and "photon plus jet", respectively. The P^"^* and 6'-''^* distributions for the 
hadronic jet with the largest transverse momentum are shown in figure[3l Both distributions are 
reasonably well described by the sum of the scaled signal and background MCs. 

A total of 14670 events with a scattered electron and an isolated photon candidate are se- 
lected, of which 6495 have at least one additional hadronic jet. 



5 Photon Signal Extraction 
5.1 Shower shape analysis 

The extraction of the photon signal exploits the fine granularity of the electromagnetic part of 
the LAr calorimeter. In order to discriminate between signal photons and the background from 
neutral hadrons and their decay products, the calorimeter cluster corresponding to the isolated 
photon candidate is further analysed using the following six shower shape variables calculated 
from the measurements of the individual cells composing the cluster [36]: 

1. The fraction of the energy of the electromagnetic cluster contained in the cell with the 
largest energy deposit ("hottest cell"). 

2. The fraction of the energy of the electromagnetic cluster contained in four or eight (de- 
pending on the granularity of the calorimeter) contiguous cells in the first two calorimeter 
layers. The cells include the hottest cell and are chosen to maximise the energy which 
they contain ("hot core"). This and the first variable are sensitive to the compactness 
of the cluster in the calorimeter. The values of these variables are on average larger for 
photons than for the background. 

3. The fraction of the cluster's energy detected in the first calorimeter layer ("layer 1"), 
which is expected to be larger on average for multi-photon clusters than for those initiated 
by a single photon. 

4. The transverse! symmetry St of a cluster defined as the ratio of the spread (defined by 
root mean squared) of the transverse cell distributions along the two principal axes. A 
photon cluster is expected to be symmetric with St values close to unity, whereas multi- 
photon clusters are typically more asymmetric and yield lower St values. 



'^In the context of the cluster shape analysis the transverse plane is defined as perpendicular to the direction of 
the photon candidate. 



5. The transverse radius of the cluster defined as the square root of the second central trans- 
verse moment R^ = ^JJ^i, where the fc'th central transverse moment of the cells distribu- 
tion is given by /i^ = (|r — {f)\^). Here r is the transverse projection of a cell position 
and (f) = (XliLo ^-'ifi)! S ^i the energy weighted average of the cell positions r^ in the 
plane transverse to the photon direction. As explained in section 4.2, only events with a 
cluster candidate of small transverse radius i?^ < 6 cm are selected for the multivariate 
analysis. 

6. The transverse kurtosis Kt is defined as Kt = /X4/(/i2)^ — 3. It specifies how strongly 
the energy distribution is peaked and is equal to zero for a Gaussian distribution. 

The discrimination power of signal and background becomes weaker at high transverse ener- 
gies, where the multi-photon clusters become more similar to a single photon cluster. Therefore 
events with E]. > 10 GeV are excluded from the measurement, as described in section 4.2. 

The distributions of the six shower shape variables are shown for the isolated photon candi- 
dates in figure lU The data are compared with the sum of the background and the scaled signal 
MC distributions. A good agreement is observed. 



5.2 Signal extraction 

In order to discriminate between single photons (signal) and single neutral hadrons (back- 
ground), probability density functions p*"^ ' are determined for the six shower shape variables, 
using simulated "single particles events", described in section 3. The signal probability density 
functions pl^ are simulated using single photon events whereas the background probability den- 
sity functions pi are approximated using a sample of events containing single neutral hadrons 
(7r°, 7], 7]', p, Lu, K*, K1, Kg, n and n). The relative contributions of the various neutral hadrons 
species are taken as predicted by the RAPGAP generator. In particular, 7r° and r] mesons con- 
tribute to 90% of the background. 

An overlap of clusters of different particles can occur due to large multiplicities specific to 
the hadronic environment in DIS. The overlap affects the photon candidate cluster shape. It is 
found to be important only for the background and leads to a loss of photon candidates due to 
the distortion of the transverse cluster radius, which then exceeds in most cases the upper limit 
of 6 cm required in the analysis (cf. section 4.2). This loss due to cluster overlap is modelled 
by supressing the background single particles according to a probability p^o proportional to the 
transverse area of the cluster pco = a--Ry witha = 0.004 cm~^. The constant a is determined by 
a comparison of single particles event samples with full MC simulation in phase space regions 
where sufficient statistics are available. 

The multi-dimensional photon and background probability densities are taken as the product 
of the respective shower shape densities P-y^bg = n«=i &P\ by ^^^ Q^ch. event a discriminator 
{D) is formed, which is defined as the photon probability density divided by the sum of the 
probability densities for photons and background D = Py/{Py + Phg)- The discriminator has 
in general larger values for isolated photons than for the decay photons. Figure [5] shows the 
discriminator distribution for the data together with the predictions of the background and the 
scaled signal MC. The data are well described by the sum of the MC predictions. 
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Since the shower shape densities vary significantly as a function of the cluster energy and 
depend on the granularity of the LAr calorimeter, the discriminator is determined in bins of 
{Erj., 7]"'), with three intervals in E], and five in rj'^. The intervals in rj"' correspond to the 
different wheels of the calorimeter. The contributions of photons and neutral hadrons in any 
of the 15 analysis bins is determined by independent minimum-^^ fits to the data discriminator 
distributions. In each (E^, r]"') bin, the x^ function is defined to be 

2 _ V^ {J^data,i ^ J^bgd-bg^ ~ J^sigdsig,i) 
j data,i ' bg bg,i ' sig sig,i 

where the sum runs over the bins of the discriminator distributions. Ndata,i is the number of data 
events in the i'th bin. dsig^i and dbg,i denote the i'th bin content of the signal and background 
discriminator distribution, respectively, normalised to unity. The ai represent the associated sta- 
tistical errors. Ngig and Nbg, represent the number of signal and background events respectively, 
and are determined by the fit. If the content in any data histogram bin is small {Ndata,i < 7), 
adjacent bins are merged. The fit is performed separately for the three samples selected as 
presented in section 4.2: inclusive, photon plus no-jets and photon plus jet. The fit quality is 
acceptable in all differential bins. 

The total number of isolated photons is obtained by summing the Nsig from all analysis 
bins. As a result 4372 ± 145 signal events are attributed to the inclusive data sample. In the 
photon plus no-jets and photon plus jet subsamples 1755 ± 106 and 2606 ± 95 signal events are 
found, respectively. The signal fraction is 29% in the inclusive sample, 20% for the photon plus 
no-jets sample and 39% for the photon plus jet sample. 



5.3 Cross section determination 

In each bin of the kinematic variables, the cross section a is computed from the number of 
events with photons in the corresponding bin as o- = Nsig/{C ■ e), where C is the luminosity. 
The correction factor e takes into account the acceptance, trigger and reconstruction efficien- 
cies, and migration between the bins. It is calculated using the scaled signal MC. Bin averaged 
cross sections are quoted in all tables and figures. The total inclusive cross section is obtained 
by summing the measured cross sections from all 15 analysis bins. The bin averaged single 
differential cross sections da/dE^ and da/drf' are obtained accordingly by summing all cor- 
responding bins in rf^ and E"^, respectively. The single differential cross section da /dQ"^ is 
determined by fitting the discriminator distributions separately in five different bins in Q^. In 
these fits, the signal and background discriminator distributions in each E^. and rj^ bin are as- 
sumed to be independent of Q'^. It has been verified that the variation of the Q"^ dependence 
have indeed a negligible effect. 



6 Systematic Uncertainties 

Two additional event samples are used for the determination of systematic errors and in-situ 
energy calibration. The first sample, containing Bethe Heitler events, ep -^ e^p, consists of 
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events with an electron reconstructed in the LAr calorimeter, a photon in the SpaCal and nothing 
else in the detector. The second, complementary, sample is selected by requiring an electron 
in the SpaCal and a photon in the LAr calorimeter and no other particle detected in the event. 
Such events originate to a large part from deeply virtual Compton scattering ep -^ e'jp. These 
independent event selections, denoted BH and DVCS respectively, provide a clean sample of 
electromagnetic clusters in the LAr calorimeter. 

The effects of the different systematic errors on the cross section are evaluated by applying 
variations to the MC simulation. The following uncertainties are considered: 

• The measured shower shape variables in the BH and DVCS event samples are compared 
to simulated single particle photons and electrons, respectively. The uncertainty on the 
shower shape simulation is estimated by distorting the discriminating variables within the 
limits deduced from the differences observed for the control samples (BH and DVCS) 
between data and simulation. The fits for the signal extraction (section 5.2) are repeated 
with the distorted distributions of the discriminating variables. The resulting systematic 
error on the total inclusive cross section is +10.2% and —12.8%. It varies between 11% 
and 25% for the single differential cross sections. The error increases with increasing E^ 
and towards large r]"' and is independent of Q^. 

• The uncertainty on the photon energy measurement is estimated using the BH and DVCS 
control samples. For the BH events the cluster transverse energy is compared with the 
track measurement. For DVCS events the energy of the cluster is compared to the energy 
calculated using the double angle method [37]. The photon energy scale uncertainty 
estimated with this method varies from 1% for photons detected in the backward region 
to 4% for forward photons. The resulting error on the total inclusive cross section is 
+3.6% and -2.6%. 

• An uncertainty of 3 mrad (4 mrad for 77-,, > 1.4) is attributed to the measurement of 
the polar angle of the photon. This uncertainty is determined by comparing the polar 
angle measurements of the track and the cluster for identified electrons in the BH control 
sample. The resulting error on the total inclusive cross section is +0.1% and —0.7%. 

• An uncertainty of up to 2% is attributed to the energy of the scattered electron and an 
uncertainty of 2 mrad to the measurement of the scattering angle. The resulting error on 
the total inclusive cross section is +1.9% and —2.9%. 

• A 3% uncertainty is attributed to the energy of hadronic final state objects [38]. This 
affects the total cross section by +1.2% and —0.7%. 

• A 5% uncertainty is applied for the correction of the cluster overlaps in the signal extrac- 
tion procedure, corresponding to half the size of the correction. 

• The uncertainty attributed to the model dependence of the acceptance corrections is de- 
rived from the differences between the acceptances determined with the PYTHIA and 
HERWIG simulations. A resulting error of 5% is found for the inclusive and the photon 
plus jet and 10% for the photon plus no-jets cross sections. 
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• The uncertainty of the ratio of the LL and QQ contributions for the acceptance correc- 
tions is taken into account by varying the scaling factor for PYTHIA from 1.5 to 3. The 
resuhing error on the total inclusive cross section is +1.4% and —0.8%. For the dou- 
ble differential cross section, a systematic error of up to 5% is found. In regions with 
7]'^ > —0.6, the systematic error is below 1%. 

• An uncertainty of 1% is attributed to the simulation of the trigger efficiency. 

• The uncertainty on the track reconstruction efficiency results in an error of ±0.3% for the 
total inclusive cross section. 

• The uncertainty on the conversion probability of the photons before entering the calorime- 
ter results in a systematic error of 2% of the cross sections. 

• The uncertainty on the luminosity measurement is 3.4%. 

In each analysis bin the individual effects of these experimental uncertainties are combined in 
quadrature. The systematic uncertainty obtained on the total inclusive cross section is +13.6% 
and —15.5%. The largest contribution to this uncertainty arises from the systematic effects 
attributed to the description of the shower shapes, which is partially correlated among measure- 
ment bins. 



7 Theoretical Predictions 



The measured cross sections are compared to the unsealed MC prediction as discussed in sec- 
tion 3. In addition, data are compared with fixed order QCD calculations, described in this 
section. A LO 0{a^a^) calculation [11,12] is used. The photon plus jet cross sections are fur- 
ther compared to a NLO 0{a^as) QCD calculation [15] which is only available for the photon 
plus jets phase space selection. In the calculations, the cross section a{ep — »■ e'jX) is obtained 
by convoluting the parton-level cross sections (for instance a{eq — > e'jq) at LO) with the proton 
parton density functions. 

The QQ contribution is dominated by the direct radiation of the photon from the quark in- 
volved in the parton level process, but also contains the contribution from quark fragmentation 
to a photon [39,40]. The direct part can be calculated in perturbation theory. The fragmentation 
contribution is described by a DIS matrix element eg -^ eq convoluted with a process indepen- 
dent quark-to-photon fragmentation function determined from data and denoted by Dq^^{z). 
Here z is the fraction of the quark's momentum carried by the photon. Already at leading order, 
the parton-level cross section contains an infrared divergence due to the emission of a photon 
coUinear to the quark. It is factorised into the fragmentation function at a factorisation scale 
jip- This singularity at LO makes a NLO prediction for the inclusive cross section considerably 
involved. 

So far, only two measurements exist that give direct information on the quark-to-photon 
fragmentation function. These measurements were made by the LEP experiments ALEPH [34] 
and OPAL [41]. Only the ALEPH measurement has the precision to allow a determination of 
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Dq^j(z). The LO calculation used in this analysis is based on the ALEPH LO parametrisation 
of D(^^^{z) [33, 34]. The NLO calculation uses the fragmentation function from [42]. In the 
MC models the contribution from fragmentation is not included. 

For photons plus no-jets the quark-to-photon fragmentation enters already at lowest order 
(a^a°) [12], in contrast to the photon plus jet sample where the fragmentation contribution is of 
the order a^a\. Since the contribution from fragmentation is suppressed by the requirement of 
the photon being isolated, the present analysis has no sensitivity for a determination oi Dq^^{z). 

In the LL subprocess, the photon is radiated by the lepton. In the present analysis, the 
reconstruction of the photon and the electron in different parts of the detector ensures that 
the two particles are separated, hence the LL term contains no coUinear singularities. The 
interference term {LQ), which differs in sign for e^p and e^p scattering, is included in the 
calculations. It contributes less than 3% to the cross section [11]. 

The calculations of the isolated photon cross section are made for the same ratio of lumi- 
nosities for e^p (47.9%) and e~p (52.1%) scattering as for the data. The cuts in the theoretical 
calculations are adjusted to correspond to the experimental cuts described in this paper. The 
calculations use the same jet algorithm, applied on parton level, as used for the data analysis. 
The CTEQ6L [29] leading order parametrisation of proton parton distributions is used. Differ- 
ent proton parton density functions are found to change the predictions by 5 — 10%. The LO 
and NLO predictions are compared to the data after a correction for hadronisation effects. The 
correction factors fhad are defined as the ratio of the cross sections calculated from hadrons to 
those from partons and are determined from the scaled signal MC. The same jet algorithm as 
for the data is applied on parton and on hadron level. The uncertainty of the correction fac- 
tor is estimated by comparing the correction factors obtained from PYTHIA with those from 
HERWIG. The correction for the total inclusive cross section is —14% with an associated un- 
certainty of 5%. For the differential cross sections the corrections are at most —30%. The 
correction is largest for low photon energies and in forward direction of the photon because of 
the degraded resolution of the isolation parameter z. The combined uncertainties of the theo- 
retical predictions from hadronisation corrections and proton parton distributions amount to up 
to 11%. 

The NLO calculation of the photon plus jet cross sections includes processes with an addi- 
tional gluon, either as the incoming parton or in the final state, as well as virtual corrections. 



The renormalisation and factorisation scales are set to /Uij= /iF = \J Q"^ + {P^'^Y- Theoretical 
uncertainties are estimated by varying jir and jip independently by a factor two up and down. 
These uncertainties are below 3% and lower than the uncertainties from the choice of the proton 
parton distributions and from the hadronisation corrections. 

8 Results 

The isolated photon cross sections presented below are given for the phase space defined in 
table ill Bin averaged differential cross sections are presented in figures [6] to [10] and in tables |2] 
to [51 For all measurements the total uncertainty is dominated by systematics. The results are 
compared with the signal MC predictions (unsealed PYTHIA plus RAPGAP) and with the LO 
and NLO calculations. The factors fhad for the correction of the theoretical calculations from 
parton to hadron level are given in the cross section tables. 
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8.1 Inclusive isolated photon cross sections 

The measured inclusive isolated photon cross section is 

a{ep -^ e-fX) = 50.3 ± 1.7 (stat) tji (syst) pb. 

The LO calculation predicts a cross section of 28.6 pb, while the signal MC expectation is 
26.4 pb. Thus both predictions significantly underestimate the measured total inclusive cross 
section by almost a factor of two. 

Differential cross sections da/dE]., da /drf' and da jdQ'^ are presented in table |2] and in 
figure [6l 

The data are compared in the left panels of figure [6] with the LO predictions, displaying 
separately the LL and QQ contributions. The LO calculation provides a reasonable description 
of the shapes of the data distributions in _E^ and rp , while the global normalisation is about a 
factor of two too low. The calculation is closest to the data at larger Q'^ and for backward 
photons (rp < —0.6), where the LL contribution is of similar magnitude to that of QQ. For 
forward and central photons (r]'^ > —0.6), the QQ contribution dominates. 

The measurements are compared in the right panels of figure [6] with the signal MC pre- 
dictions. The estimations of the QQ processes by PYTHIA and of the LL contributions by 
RAPGAP agree well with the predictions from the LO calculation. Thus the conclusions for the 
comparison of the signal MC with the data are similar as for the LO calculation. This agreement 
between LO and the signal MC holds for all cross sections studied in this analysis and in the 
following the data are only compared with the LO calculation. 

Figure|7]and table [3] present the differential cross section as a function of Ej, in five different 
bins of r]"^, corresponding to the wheel structure of the LAr calorimeter. The shapes of all cross 
section distributions are reasonably well described by the LO calculation. A good description of 
the data can be obtained by a global scaling of the QQ contribution by a factor 2.3 and leaving 
the LL normalisation unchanged. This is an indication that the observed excess in the data is 
mainly due to an underestimation of the QQ component. 

At Q^ > 40 GeV^, the agreement of the predictions with the data is somewhat better, as can 
be seen in figure [8] and table lU The total inclusive cross section for Q'^ > 40 GeV^ is 

a{ep -^ e-fX) = 14.0 ± 0.8 (stat) tH (syst) pb. 

The LO prediction of 10.3 pb is about 30% below the data. The shapes of the da/dE^ and 
da/dr]"' distributions are well reproduced. The relative contribution of LL is predicted to be 
higher than at low Q^. 

The present measurement is extrapolated to the phase space of the analysis performed by 
the ZEUS collaboration [10] (Q^ > 35 QeV^, y > 0, E^ > 10 GeV, 139.8 < Oe < 171.9° 
and 5 < E^ < 10 GeV). Due to the different angular coverage of the calorimeters of the 
two detectors, the extrapolation is not possible for r]^ < —0.6 because the photon-electron 
separation cannot be properly taken into account. Figure |9] shows the comparison of the HI and 
ZEUS mesurements of the differential cross section da/dr]'^. A good agreement is observed. 
The LO calculation is also shown and is here corrected for hadronisation effects in contrast 
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to the comparison in [11]. In the ZEUS analysis the photon radiation from the electron is 
neglected in the acceptance corrections and no Wx cut is applied. A rough estimate shows that 
with the Wx cut used by HI and the acceptance calculation with the combination of PYTHIA 
and RAPGAP, the ZEUS cross section values are expected to be lowered by about 10 — 30%. 
A different photon isolation criterion is used in the ZEUS analysis which is not corrected for. 
Studies of isolated photons in photoproduction indicate that the two different isolation criteria 
give very similar results. The measurement presented in this paper significantly extends the 
kinematical region probed by ZEUS in photon transverse energy and pseudorapidity, and in Q^. 



8.2 Photon plus no-jets and photon plus jet cross sections 

The cross section for jet production in events with isolated photons is studied. The measurement 
is performed in the phase space defined for the inclusive cross section with an additional jet 
requirement or veto as shown in table [B 

The measured total cross section for photon plus no-jets is 

a{ep -^ e-fY) = 18.8 ± 1.2 (stat) tH (syst) pb, 

where Y contains no identified hadronic jet. The prediction of the LO calculation is 11.7 pb. 
The measured total cross section for the photon plus at least one jet is 

a{ep -^ e-fjet X) = 31.6 ± 1.2 (stat) ItsCsyst) pb. 

The prediction of the LO calculation is 16.7 pb. For both samples the predictions are sig- 
nificantly lower than the data. The observed ratios of data to LO prediction are very similar to 
the inclusive isolated photon measurement. As for the inclusive sample, similar conclusions are 
found for the MC predictions. 

A comparison to a NLO calculation is possible for the photon plus jet cross sections. The 
NLO calculation predicts a cross section of 20.2 ± 0.6 pb, about 20% higher than the LO 
prediction but still roughly 35% lower than the data. The analysis performed using a higher cut 
on P^^ > 4 GeV lead to a similar disagreement between the data and the calculation. 

The measured differential cross sections for the photon plus no-jets and photon plus jet 
selections are presented in table [5l They are compared with the LO predictions in figure [lOl For 
both samples the LO prediction describes the shapes of the da/dE^ and da/drj'^ distributions 
reasonably well and is furthest below the data at low Q"^, where the QQ term dominates. All 
four diagrams in figure [U contribute to the photon plus no-jets and photon plus jet samples, but 
the LL contribution is largely suppressed for the photon plus no-jets sample due to the cut on 
Wx- Since, at leading order 0(a°), the quark fragmentation contribution to the cross section 
enters only the photon plus no-jets sample, the observed excess can not solely be attributed 
to an underestimation of that contribution. The cross section for photon plus jet production 
is roughly two times higher than for photon plus no-jets. This is in contrast to the inclusive 
ep — > eX cross section, where topologies with an additional jet are suppressed by 0{as). The 
similar cross sections for photon events with or without additional jets can be explained by the 
fact that both topologies correspond to the same order 0{a^) in perturbative QCD. 
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In addition, the differential cross sections da/dE'^, da/dr]"' and da/dQ'^ for the photon 
plus jet selection are compared to the NLO prediction (figure [10] right). On average, the NLO 
prediction is higher than the LO prediction, in particular at low Q'^, but is still lower than the 
data by roughly 35%. The shapes of all three differential cross sections are described well by 
the NLO prediction. 



9 Conclusions 

The cross section for isolated photon production ep — > e'jX is measured in deep-inelastic 
scattering at HERA. Photons with a transverse energy in the range 3 < E^ < 10 GeV and 
with pseudorapidity —1.2 < r]"^ < 1.8 are measured in DIS events in the kinematic region 
4: < Q^ < 150 GeV^, y > 0.05 and Wx > 50 GeV. Compared to previous measurements, 
the range of Q'^ is extended from Q^ > 35 GeV^ to Q^ > A GeV^. The cross section receives 
contributions from photon radiation by the struck quark (QQ), as well as from wide angle 
bremsstrahlung of the initial and final electrons (LL) and their interference. 

The data are compared to a LO 0{a^a^) calculation which is shown to underestimate the 
measured cross sections by roughly a factor of two. The prediction is most significantly below 
the data at low Q^. The shapes of the da/dE^ and da/drf^ distributions are described reason- 
ably well. The comparison of data to the LO calculation in bins of t]'^ show that the difference 
in normalisation can mainly be attributed to an underestimation of the QQ contribution. The 
data are further compared to predictions from the MC generators PYTHIA, for the simulation of 
photons radiated from the quark, and RAPGAP for photons radiated from the electron. The MC 
predictions are very similar to the predictions from the LO calculation and also understimate 
the data, especially at low Q^. 

Jet production in events with isolated photons is also studied. The cross sections for events 
with no or at least one hadronic jet are underestimated by the LO prediction by a similar factor 
as for the inclusive measurement. Again the expectations are furthest below the data at low Q'^. 
The total photon plus jet cross section is roughly double the photon plus no-jets cross section 
as expected from the calculations. The NLO 0{a^a\) prediction for photon plus jet is higher 
than the LO prediction, most significantly at low Q"^, but still underestimates the data. The NLO 
calculation describes the shapes of the differential cross sections reasonably well. 

Further theoretical investigations are needed to understand the observed discrepancy be- 
tween the measurements and the predictions, including for instance the calculation of higher 
order processes. 
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Isolated Photon Cross Section Phase Space 


Inclusive 
cross section 


3 < EJ < 10 GeV 
-1.2 <r/T < 1.8 

z = E^/E^^°*°"-^'^* > 0.9 
Ee > 10 GeV 
153 < ^e < 177° 

4 < g2 < 150 GeV2 

Wx > 50 GeV 

y > 0.05 

y^ = 319 GeV 


Jet definition 


kr algorithm with Py-weighted 
recombination scheme [35], -Ro = 1 
Pf'* > 2.5 GeV 

-1.0 < r/J"* < 2.1 (hadronicjet) 
-2.0 < r^p^oto'^-ie* < 2.1 (photon-jet) 



Table 1: Phase space region in which isolated prompt photon cross sections are measured 
together with the definition of jets. Kinematics are defined in the HI laboratory frame. 
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HI Inclusive Isolated Photon Cross Sections 


EJ 


da/dEl 


Stat. 


syst. 


Jhad 


[GeV] 


[pb/GeV] 






3.0 4.0 


16.98 


±1.20 


+2.79 
-2.61 


0.78 


4.0 6.0 


10.51 


±0.47 


+1.50 
-1.86 


0.89 


6.0 10.0 


3.08 


±0.20 


+0.46 
-0.60 


0.98 


^7 


da/drj'^ 


Stat, 
[pb] 


syst. 


fhad 


-1.2 -0.6 


26.15 


±1.67 


+3.60 
-4.16 


0.92 


-0.6 0.2 


20.69 


±1.34 


+3.53 
-3.73 


0.85 


0.2 0.9 


15.83 


±0.93 


+1.97 
-3.25 


0.81 


0.9 1.4 


9.57 


±0.87 


+1.99 
-2.00 


0.80 


1.4 1.8 


5.50 


±1.15 


+1.04 
-1.75 


0.80 


Q' 


da/dQ^ 


Stat. 


syst. 


Jhad 


[GeV2] 


[pb/GeV2] 






4.0 10.0 


2.48 


±0.21 


+0.34 
-0.41 


0.87 


10.0 20.0 


1.17 


±0.07 


+0.19 
-0.21 


0.83 


20.0 40.0 


0.52 


±0.03 


+0.07 
-0.10 


0.81 


40.0 80.0 


0.235 


±0.013 


+0.033 
-0.048 


0.83 


80.0 150. 


0.063 


±0.006 


+0.009 
-0.012 


0.87 



Table 2: Differential cross sections for inclusive isolated photon production in the kinematic 
range specified in table [IJ fhad denotes the hadronisation correction factor applied to the LO 
calculation. 
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HI Inclusive Isolated Photon Cross Sections 


h/rp 


da/dEl 


Stat. 


syst. 


J had 


[GeV] 


[pb/GeV] 










-1.2 < 


r/^ < -0.6 






3.0 


4.0 


4.86 


±0.67 


+0.88 
-0.63 


0.86 


4.0 


6.0 


3.46 


±0.28 


+0.48 
-0.66 


0.96 


6.0 


10.0 


0.98 


±0.12 


+0.13 
-0.23 


1.00 






-0.6 < 


: T]"' < 0.2 






3.0 


4.0 


5.81 


±0.75 


+1.27 
-1.20 


0.76 


4.0 


6.0 


3.20 


±0.28 


+0.56 
-0.65 


0.88 


6.0 


10.0 


1.09 


±0.13 


+0.15 
-0.17 


0.99 






0.2 < 


rf < 0.9 






3.0 


4.0 


3.94 


±0.51 


+0.59 
-0.72 


0.72 


4.0 


6.0 


2.39 


±0.16 


+0.28 
-0.51 


0.84 


6.0 


10.0 


0.59 


±0.06 


+0.09 
-0.14 


0.96 






0.9 < 


T]^ < 1.4 






3.0 


4.0 


1.66 


±0.31 


+0.22 
-0.31 


0.69 


4.0 


6.0 


0.82 


±0.12 


+0.21 
-0.16 


0.82 


6.0 


10.0 


0.37 


±0.05 


+0.10 
-0.09 


0.96 






1.4 < 


r]^ < 1.8 






3.0 


4.0 


0.72 


±0.28 


+0.23 
-0.21 


0.70 


4.0 


6.0 


0.64 


±0.16 


+0.09 
-0.22 


0.81 


6.0 


10.0 


0.049 


±0.045 


+0.009 
-0.016 


0.94 



Table 3: Differential cross sections for inclusive isolated photon production da/dE^ in differ- 
ent 7]'^ bins corresponding to the wheel structure of the LAr calorimeter (see text). The kinematic 
range is defined in table \T\ fhad denotes the hadronisation correction factor applied to the LO 
calculation. 
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HI Inclusive Isolated Photon Cross Sections 


for Q2 > 40 GeV^ 






h/q^ 


da/dE^ Stat. 


syst. 


fhad 


[GeV] 


[pb/GeV] 






3.0 4.0 


3.70 ±0.39 


+0.63 
-0.59 


0.80 


4.0 6.0 


2.53 ±0.23 


+0.43 
-0.43 


0.87 


6.0 10.0 


1.30 ±0.15 


+0.19 
-0.27 


0.96 


^7 


da/dif Stat. 
[pb] 


syst. 


Jhad 


-1.2 -0.6 


9.61 ±1.00 


+1.48 
-2.02 


0.97 


-0.6 0.2 


5.13 ±0.59 


+0.92 
-0.97 


0.86 


0.2 0.9 


3.49 ±0.32 


+0.45 
-0.72 


0.78 


0.9 1.4 


2.37 ±0.33 


+0.49 
-0.47 


0.76 


1.4 1.8 


1.12 ±0.51 


+0.21 
-0.34 


0.73 



Table 4: Differential cross sections for inclusive isolated photon production in the kinematic 
range specified in table [Hand 40 < Q^ < 150 GeV^. fhad denotes the hadronisation correction 
factor applied to the LO calculation. 
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HI Photon plus no- Jets 


HI Photon plus Jet 


E^ 


da/dE:^ 


Stat. 


syst. 


fhad 


da/dE:{, 


Stat. 


syst. 


fhad 


[GeV] 


[pb/GeV] 






[pb/GeV] 






3.0 4.0 


8.10 


±0.93 


+1.82 
-1.48 


0.75 


8.85 


±0.70 


+1.37 
-1.38 


0.82 


4.0 6.0 


3.79 


±0.29 


+0.69 
-0.77 


0.91 


6.65 


±0.35 


+0.92 
-1.15 


0.89 


6.0 10.0 


0.77 


±0.10 


+0.14 
-0.18 


1.10 


2.35 


±0.17 


+0.35 
-0.46 


0.97 


rj-y 


da/drj^ 


Stat, 
[pb] 


syst. 


Jhad 


da /dr]'^ 


Stat, 
[pb] 


syst. 


fhad 


-1.2 -0.6 


9.30 


±1.07 


+1.53 

-1.72 


0.88 


16.61 


±1.20 


+2.67 
-2.64 


0.97 


-0.6 0.2 


8.46 


±0.95 


+1.73 
-1.84 


0.81 


12.32 


±0.90 


+2.15 
-2.19 


0.88 


0.2 0.9 


5.98 


±0.71 


+1.16 
-1.34 


0.82 


9.94 


±0.59 


+1.16 
-2.03 


0.81 


0.9 1.4 


2.57 


±0.47 


+0.64 
-0.61 


0.85 


6.99 


±0.73 


+1.38 
-1.40 


0.79 


1.4 1.8 


2.40 


±0.73 


+0.78 
-0.74 


0.91 


3.22 


±0.85 


+0.61 
-1.01 


0.77 


Q' 


da/dQ^ 


Stat. 


syst. 


j had 


da/dQ^ 


Stat. 


syst. 


fhad 


[GeV2] 


[pb/GeV2] 






[pb/GeV2] 






4.0 10.0 


1.09 


±0.16 


+0.21 
-0.22 


0.88 


1.39 


±0.13 


+0.20 
-0.22 


0.87 


10.0 20.0 


0.44 


±0.05 


+0.09 
-0.09 


0.81 


0.76 


±0.05 


+0.11 
-0.14 


0.86 


20.0 40.0 


0.21 


±0.02 


+0.04 
-0.05 


0.80 


0.31 


±0.02 


+0.04 
-0.06 


0.83 


40.0 80.0 


0.071 


±0.008 


+0.012 
-0.017 


0.81 


0.162 


±0.010 


+0.024 
-0.033 


0.84 


80.0 150. 


0.021 


±0.004 


+0.005 
-0.007 


0.88 


0.040 


±0.005 


+0.005 
-0.007 


0.89 



Table 5: Differential cross sections for the production of isolated photons accompanied by 
no or at least one hadronic jet in the kinematic range specified in table [T] f^ad denotes the 
hadronisation correction factor applied to the LO and the NLO calculation. 
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Figure 2: Distributions of (a) iJj and (b) 6'^ for isolated photon candidates in the final event 
sample. Data are shown as points with error bars. The bold solid histogram shows the sum of the 
expectation from RAPGAP (non-rad.) for neutral hadron background (shaded), from PYTHIA 
for radiation from the quark scaled by a factor of 2.3 (dashed line) and from RAPGAP (rad.) for 
radiation from the electron (solid line). The unshaded area corresponds to the estimated isolated 
photon contribution (RAPGAP(rad) plus PYTHIA x 2.3). 
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Figure 3: Distributions of (a) the transverse momentum and (b) the polar angle of the hadronic 
jet with the highest transverse momentum in events with an isolated photon candidate. The data 
are compared to the MC predictions, as described in the caption of figure |2l 
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Figure 4: Distributions of the six variables that are used to define the discriminant for isolated 
photon identification: (a) hottest cell fraction, (b) fraction of the hot core, (c) first layer fraction, 
(d) transverse symmetry St, (e) transverse radius Rt and (f) transverse kurtosis Kt- The data 
are shown with the MC predictions described in the caption to figure [21 The shape difference 
between RAPGAP (rad.) and PYTHIA arises from the different distributions in phase space 
(see figure |2b). 
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Figure 5: The distribution of the discriminator (D) used in the identification of isolated photon 
candidates for events that have passed the event selection. The data are compared to the MC 
predictions described in the caption of figure |2l 
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H1 Isolated Photon Production in DIS 
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Figure 6: Differential cross sections for inclusive isolated photon production da/dEj,, da /drf , 
and da /dQ"^ in the kinematic range specified in table [B The inner error bars on the data points 
indicate the statistical error, the full error bars contain in addition the systematical errors added 
in quadrature. The cross sections in (a, c, e) are shown together with a leading order, C(a^a°), 
calculation corrected for hadronisation effects, LL corresponding to radiation from the electron, 
QQ to radiation from the quark and LQ to their interference. The same cross sections are 
shown in (b, d, f) together with the prediction from PYTHIA for photon emission from the 
quark and from RAPGAP (rad.) for emission from the electron. 
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H1 Isolated Photon Production in DIS 
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Figure 7: Differential cross sections da/dE]. for inclusive isolated photon production in the 
kinematic range specified in table [U in rj"' bins corresponding to the wheel structure of the LAr 
calorimeter (see text). The measurements are compared to a leading order C(a^a°) calculation 
(more details in the caption of figure[6l). The LL contribution is negligible for r]^ > 0.2. 
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H1 Isolated Photon Production in DIS at Q^ > 40 GeV^ 
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Figure 8: Differential inclusive cross sections da/dEj, (a) and da /drf' (b) for isolated photon 
production in the kinematic range specified in table [Tj and the additional criterion 40 < Q^ < 
150 GeV^. The cross sections are compared to the leading order, C(a^Q;°), calculation (more 
details in the caption of figure [6]). 
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Figure 9: Differential cross sections da /drf^ for the inclusive isolated photon production in 
comparison to the previous measurement by ZEUS [10] for Q"^ > 35 GeV^, E'^ > 10 GeV, 
139.8 <e^< 171.9° and 5 < Ej < 10 GeV. The additional condition Wx > 50 GeV is used 
in the HI analysis only (see section 8). The cross sections are compared to the leading order, 
0{a^a^), calculation. In contrast to the comparison in [11], the calculation is here corrected 
for hadronisation effects. 
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Figure 10: Differential cross sections da /dE^, da /drf' and da/dQ"^ for photon plus no-jets (a, 
c, e), and photon plus jet (b, d, f) production in the kinematic range specified in table [IJ The 
cross sections are compared to the leading order, C(a'^a°), calculation (more details in the 
caption of figure [6l). The photon plus jet sample is additionally compared to a NLO {a^ag) 
calculation. The bin averaged NLO cross sections are indicated by the squares. 
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